Theory of ground state cooling of a mechanical oscillator using dynamical back-action 
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A quantum theory of cooling of a mechanical oscillator by radiation pressure-induced dynamical 
back-action is developed, which is analogous to sideband cooling of trapped ions. We find that final 
occupancies well below unity can be attained when the mechanical oscillation frequency is larger 
than the cavity linewidth. It is shown that the final average occupancy can be retrieved directly 
from the optical output spectrum. 
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Mesoscopic mechanical oscillators are currently at- 
tracting interest due to their potential to enhance the 
sensitivity of displacement measurements [l| and to 
probe the quantum to classical transition of a macro- 
scopic degree of freedom @, EJ . A prerequisite for these 
applications is the capability of initializing an oscilla- 
tor with a long phonon lifetime in its quantum ground 
state. So far this has not been demonstrated because 
the combination of sufficiently high mechanical frequen- 
cies (w m /27r) and quality factors in the relevant regime 
hw m 3> ksT has not been reached In contrast, in 
atomic physics laser cooling has enabled the preparation 
of motional ground states [4|, l5( . This has prompted re- 
searchers to study means of cooling a single mechani- 
cal resonator mode directly using laser radiation. Early 
work demonstrated cooling of a mechanical degree of free- 
dom of a Fabry-Perot mirror using a radiation pressure 
force controlled by an electronic feedback scheme @, 0] , 
in analogy to stochastic cooling. In contrast, the radia- 
tion pressure induced coupling of an optical cavity mode 
to a mechanical oscillator [cf. Fig. Ufa)] can give rise to 
self-cooling via dynamical back-action Q. In essence, 
the cavity delay induces correlations between the radi- 
ation pressure force and the thermal Brownian motion 
that lead to cooling or amplification, depending on the 
laser detuning. In a series of recent experiments, these 
effects have been used to cool a single mechanical mode 
S EH E3 • While classical and semiclassical analysis of 
dynamical back-action have been developed ID, 14|, the 
question as to whether ground state cooling is possible 
has not been addressed. 

Here a quantum theory of cooling via dynamical back- 
action is presented. We find that final occupancies below 
unity can indeed be attained when the optical cavity's 
lifetime is comparable to or exceeds the mechanical os- 
cillation period. Along these lines, an analogy between 
this mechanism and the sideband cooling of trapped ions 
in the Lamb-Dicke regime is elucidated [5| . In our setting 
the optical cavity mode plays the role of the ion's pseu- 
dospin mediating the frequency up-conversion underlying 
the cooling cycle. Finally, we discuss how the average 
phonon occupancy can be retrieved from the spectrum 
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FIG. 1: (color online) (a) Fabry-Perot equivalent of a mechan- 
ical eigenmode (frequency uj m /2n and quality factor Q m ) cou- 
pled to an optical cavity mode, (b) Processes that decrease 
or increase the mechanical eigenmode's quantum number n. 
Cooling occurs by scattering into anti-Stokes modes, whereas 
heating (amplification) proceeds via Stokes scattering, (c) Il- 
lustration that compares the output power spectral density 
(solid lines), exhibiting motional sidebands (blue and red), 
with the cavity absorption in the absence of optomechanical 
coupling (dashed). 



of the optical cavity output. We note that these results 
can be applied to a wide range of experimental realiza- 
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tions of cavity self-cooling 

We treat the laser driven optical cavity mode coupled 
to the mechanical resonator mode as an open quantum 
system and adopt a rotating frame at the laser frequency 
lu l . The system Hamiltonian is given by 1(| 17 1 



H' = — HA' L a p a p + hriui m a p a p (a % 
+ (a p + a p ) + hw m a\ n a m 



(1) 



Here a p (a m ) is the annihilation operator for the op- 
tical (mechanical) oscillator, lu p (w m ) its angular fre- 
quency and A' L the laser detuning from the optical res- 
onance. We have also introduced the driving amplitude 
ft = 2y / P/haJLT cx , where P is the input laser power and 
1/t ox the photon decay rate into the associated outgo- 
ing modes (e.g. optical fiber modes [lH). The optome- 
chanical coupling via radiation pressure can be character- 
ized by the dimensionless parameter 77 = (oj p /oj m )(l m /L); 
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with l m = ^Ti/2mu> m the zero point motion of the me- 
chanical resonator mode, m its effective mass, and L an 
effective length that depends on how the radiation pres- 
sure force affects the optical cavity. For typical materials 
and dimensions 11 1 one obtains r\ ~ 10~ 4 . 

The optical cavity losses and the intrinsic dissipation of 
the mechanical resonator are characterized, respectively, 
by the cavity lifetime r and the mechanical quality fac- 
tor Q m . These give rise to a dissipative contribution to 
the Liouvillian £ D (he. £ = -{[H', . . .} + C'n) that 
is of Lindblad form [25| with collapse operators [13, EH : 
y /f T/ra p , N /7 m n(u; m )aJ n and ^"f m [n(u m ) + l]a m . Here 
7m = ^m/Qm is the mechanical oscillator's natural 
linewidth and n(u m ) its Bose number at the environmen- 
tal temperature. We will focus on the regime [n(ui m ) + 
l]7 m < uj m , j m < 1/r and ri,r)\a\ < 1, l/w m r. The 
first condition is necessary for ground state cooling (see 
below) , the second one is satisfied in all recent experi- 
ments @, H3, Uli and the last one, given the smallness of 
77, will hinge on having a sufficiently low input power. 

To study the dynamics generated by the Liouvillian 
£ it proves useful to apply a shift to the modes' nor- 
mal coordinates: a p — ► a p + a, a m — > a m + with the 
c-numbers a and (3 chosen to cancel out all the linear 
terms in the transformed Liouvillian To zeroth 

order in the small parameters r\ and 1 / Q m we have [26| : 
a VLt/(2tA' l + i), (3 w -t?H 2 . We include the ra- 
diation pressure induced optical resonance shift into the 
effective detuning A' L + 2rj 2 \a\ 2 oj m — > Al and perform 
the additional canonical transformation a p — > (a/\a\)a p . 
While the dissipative part of the Liouvillian £ remains 
invariant, the Hamiltonian transforms into 



H = - hA L a p a p + huj m a\ n a m + hriuj m [a' p a p 
+ \a\(a p + al)] (a m + a] n ) . 



(2) 



Henceforth we will refer to the primed representation 
([T|) as the "physical" one and to the unprimed represen- 
tation © as the "shifted" one. The smallness of rj 2 and 
[n(ui m ) + i\/Q m imply a wide separation between the 
timescales for cooling and heating the mechanical oscil- 
lator and those characterizing the dynamics of the opti- 
cal cavity mode and the mechanical oscillation period. 
Thus, the electromagnetic environment (including the 
optical cavity) can be regarded as a structured reservoir 
with which the mechanical mode interacts perturbatively 
[cf. Fig. QJb),(c)]. This prompts us to derive a "general- 
ized quantum optical" master equation for the reduced 
density matrix [18j of the latter: /i = Tr p {p}. In our 
context, such a derivation can also be viewed as an adia- 
batic elimination (l8j of the optical cavity in the presence 
of fast rotating terms (oc e ±luJmt ) in the optomechanical 
interaction. We note that while in the physical represen- 
tation the steady state average occupancy of the optical 
cavity is given by \a\ 2 , in the shifted one its steady state 



is simply the vacuum |0) p . Thus, we obtain 

fi = -i [w TO a m a m ,/x] + \{y m [n (u m ) + 1] + A-} 
x (2a m /m TO - a ] m a n ji - na\ n a m ) + \[^ m n (u m ) 
+ A+] (2a m ^a m - a m a^/i - ^a m o) m ) . (3) 

In the first term we have redefined u> m to include the 
light-induced shift of the mechanical frequency. The sec- 
ond and third terms correspond, respectively, to cooling 
and heating induced by the coupling to the thermal bath 
(contributions oc 7 m ) and by inelastic laser light scatter- 
ing processes [cf. Fig. Wfb)] with rates 
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In the shifted representation it is simple to understand 
these cooling and heating rates in terms of perturbation 
theory in the small parameters 7y|a| and r\ [cf. Eq. To 
lowest order in r\ only the states |0) p and |l) p participate 
yielding the same results as for an equivalent dissipative 
two level system. Denoting by |n) the number states 
of the mechanical oscillator we have anti-Stokes (Stokes) 



processes in which the transition |0) p |n) — > 
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(|0) p |n) — » 1 1) p| n+l)) followed by the decay of the cavity 
photon leads to cooling (amplification). This scenario is 
thus similar to the laser cooling of a trapped ion in the 
Lamb-Dicke regime 0, 0, [l!|, or of a nanomechanical 
resonator coupled to an "artificial atom" [2(| or an ion 
2l| . An important caveat in this analogy is that there is 
no external driving for r\ = 0. Furthermore, though the 
parameter rj 2 will play a role reminiscent of the Lamb- 
Dicke parameter — determining for example the relative 
weights of the sidebands — the efficiency of the cooling 
process will depend solely on r] 2 \a\ 2 and Eqs. (fj| will 
remain valid for arbitrary f2 provided rj 2 is sufficiently 
small. This absence of a "direct" driving amplitude also 
implies that the cubic term in Hamiltonian ([2]) does not 
contribute to the master equation (|3|) as it only generates 
terms that are higher order in rj 2 . Thus to lowest order 
there is no "diffusive channel" and the theory is equivalent 
to a quadratic Liouvillian [271 ] . 

Henceforth we focus on the regime < for which 
there is a net laser cooling rate Y = A_ — A + > 0. 
In this regime Eq. §5§ has a well defined steady state 
that transformed back to the physical representation 
yields a shifted thermal state. The corresponding steady 
state average occupancy, to which the system converges 
on the timescale l/(r + j m ), is given by (aJ Tl a m )ss = 



n f 



f \/3\ 2 with n/ 
= V 2 n 4 T 4 /(4r 
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[7 m n(w m ) + A+]/(7 m + T) and 
+ l) 2 . However the final tem- 



perature should be defined in terms of n/ as the other 
contribution arises from a coherent shift. In fact, un- 
der the conditions underpinning our approximations, the 
latter could be undone by a "slow" switch-off of the laser. 

As we start from thermal equilibrium, initially the 
number of phonons is given by 7ij = n(w m ). Thus, from 
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the expression for n t it is clear that for appreciable cool- 
ing (i.e. n/ < «i) we need T ^> 7 m . In this regime 



with 



n f 



. r 
r 



1 



4r 2 (Ax + u m y + 1 
16r 2 a, m (-A L ) 



(5) 



(6) 



Note that from Eq. ([4]) it follows that T is exactly pro- 
portional to the input power P. Since there are two 
terms in Eq. |5]), there are two regimes depending on 
which contribution dominates the behavior of nf . In the 
first regime heating is dominated by the intrinsic dissi- 
pation of the mechanical resonance and the final aver- 
age occupancy is proportional to the initial one. This 
behavior has already been demonstrated experimentally 
H 
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12j. On the other hand for sufficiently low 7 m 
and input laser power, the heating is dominated by the 
scattering of laser light. In this regime the final occu- 
pancy is given by n/ (cf. Fig. [2]). Thus the optimal 
value of n/ solely depends on the product u> m T and is 
found by minimizing with respect to the normalized de- 
tuning S = tAl (i.e. it is power independent). This 
yields <5 op t = — -y/l-j- 4a; 2 n T 2 /2, implying the fundamen- 
tal temperature limit: 

nTL = min {fif (5)} 



\(VT 



l/4a, 2 



(7) 



The regime w m r <C 1 is in essence the adiabatic limit [la ], 
since the cavity dynamics is much faster than the me- 
chanical oscillator's period. Several recent experiments 
fall into this regime 

0,113 ■ Expanding Eq. © we ob- 
tain nxL ~ l/4w m r 3> 1 precluding ground state cooling 
in this parameter regime. The corresponding final tem- 
perature is of order h/ksT in complete analogy with the 
Doppler limit of the laser cooling of harmonically bound 
atoms, where the sidebands are not resolved [H]. 

We turn now to the regime where retardation effects 
become significant (i.e. u m > 1/r). This has indeed been 
observed in recent experimental work pertaining to both 
amplification 14 1 and cooling [l(| EH of a mechanical os- 
cillator mode. In this regime the optical cavity field can- 
not respond instantaneously to the mechanical motion 
and an entirely viscous radiation pressure force may arise 
11]. In the frequency domain, this can be interpreted as 
having mechanical frequencies which exceed or are equal 
to the cavity linewidth [cf. Fig. QJc)]. Accordingly, the 
asymmetry in the Stokes and anti-Stokes scattering rates 
becomes more pronounced leading into the analog of the 
"resolved sideband" limit of the laser cooling of harmon- 
ically bound atoms. More precisely for 4w 2 „r 2 ^> 1, 
Eq. ([7]) yields tjtl ~ l/16w 2 „r 2 <C 1 implying that in this 
limit one can reach arbitrarily low temperatures (ground 
state cooling). In a concrete realization, a sound bench- 
mark to evaluate the cooling performance is whether oc- 
cupancies below unity can be attained. Eq. ([7]) leads to 
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FIG. 2: (color online) Final (steady state) average phonon 
number n/ as a function of normalized laser detuning (Alt) 
and normalized mechanical angular frequency (uj m T). The 
contour-lines indicate the values hf = 1 and 10, respectively. 
Ground state cooling is only possible in a finite detuning win- 
dow and for uj m T > l/-\/32- 



the following criterion [281 ]: titl < 1 ^mi~ > l/v32; 
and Eq. |6]) implies that within this regime occupancies 
below unity are only possible for a certain "detuning win- 
dow" - - l/4r 2 <A L + 3tu m < - l/4r 2 . 

Finally, we consider the impact of the intrinsic dissi- 
pation on the optimal value of n/ in the regime w TO r > 
1/V32 (analogous considerations can be done for the 
opposite regime). The situation is reminiscent of the 
"atomic" laser cooling of nanoresonators [13, [1H where 
the finite Q m also plays a crucial role. However, in the 
present context, the analysis is simpler given that the op- 
timal laser detuning to maximize the cooling rate T [and 
thus minimize the first term in Eq. ([5])] is of the same 
order as 5 op t ■ Hence, the only relevant issue is the upper 
bound on P required by the wide timescale separation 
underpinning our adiabatic treatment of the cooling and 
heating processes. Given the quadratic nature of our ap- 
proximate theory, already discussed, the Bose enhance- 
ment of the resonator mode plays no role and the adia- 
batic requirement reads A_ <C 1/r. Aside from the lim- 
itations of our approximate treatment, heuristic consid- 
erations imply that the timescale over which (a^a m )(t) 
reaches its steady state value can never be shorter than 
r. Thus, our treatment provides an upper bound for the 
ultimate final temperature when the finite Q m is consid- 
ered. As an illustration we consider the parameters of 
Ref. [ll| (i.e. A l t > 0.5, w m /2?r = 60 MHz, r = 3 ns). 
For a reservoir temperature of 4 K we have n,; w 1390. 
If we consider the improvements in the mechanical Q- 
values of toroid microcavities due to vacuum operation 
(Qm = 30, 000) a cooling rate of 2.8 MHz is then required 
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to reach nij m /T < 1 [cf. first term in Eq. [291 ]. 

The cooling process gives rise to photons which have 
frequencies that differ from the pump laser (lul)- Thus it 
can be studied in an experiment by measuring the spec- 
trum of the scattered light. As depicted in Fig. [TJa), we 
consider a one-sided cavity and the relevant observable is 
the output power. The input-output formalism implies 
that in the physical representation its spectrum S(u>) is 
given by the Fourier transform of e lULT {[y/l/T^ x aj ) (t+T)+ 

a L(* + T )][v /l / r exa p (0 + a in W])ss- In the shifted rep- 
resentation a p (t) — > a p (t) + a and the classical input just 
adds a c-number to the cavity steady state amplitude. 
Along the lines of our derivation of Eq. ([3]); a p (i), a p (t) 
are treated as environment operators to be reduced to 
the system operators a m (t), a m (i) by integrating out the 
corresponding Heisenberg equations of motion. A calcu- 
lation based on perturbation theory and the theory of 
quantum Markov processes [18| then yields [22| : 
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A+ (n f + 1) 
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* (lo-lu l + uj m y + Y J 

where the relative order of the corrections is given by rj 2 
for all frequencies, and we have normalized to number 
of photons per unit time and unit frequency. Here A-^ 
and rif are given respectively by Eqs. ((4]) and ([5]) and 
7 c ff = 7 m + r. As expected there is emission of blue- 
shifted (Anti-Stokes) photons associated to cooling and 
red-shifted (Stokes) photons associated to heating. These 
motional sidebands have a linewidth determined by the 
effective damping rate 7 e ff and weights (iVzp ) determined 
by the cooling and heating rates; namely, N- = —A—Uf 
and N + = ^-A + (n f + 1). 

The final occupancies can be retrieved by comparing 
the above spectra [Eq. ©] for different input powers. 
The quantity [N+(P) + N^(P)}P /[N + (P Q ) + N_(P )]P 
provides an upper bound for the ratio nj/m. Here we 
have introduced a "reference" low power Po for which 
T, A + -C 7 m implying n/(Po) ~ rii, and assumed that the 
input power P induces appreciable cooling [i. e. rif(P) -C 
rii]. It is important to note that (given n,) this upper 
bound provides an accurate direct measurement of the 
final temperature for n/ 1/2. On the other hand the 
worst case scenario occurs for Lo m T 3> 1 and nf ~ fif 
where it yields 2n/. However for nf < 1 an accurate 
measurement is afforded by the quantity 



N_(P)N + (P ) 



nj{P) 



Crii + PA^ 



with C = hw L T cx j m (4T 2 A 2 L + l)/4r 2 ?7 2 . Thus, the high 
power limit of Eq. ([9]) provides a clear signature for 
ground state cooling when it can be achieved. This is 
in stark contrast to the case of a laser cooled trapped ion 
where a well defined thermal reservoir associated to the 
intrinsic dissipation is lacking and detailed balance yields 
A r _ = N + . In a realistic scenario the central peak can- 
not be regarded as a delta function and the small relative 
weight of the sidebands poses an experimental challenge. 
This could be overcome by combining state of the art low 
noise lasers and high resolution spectroscopy with a suit- 
able lock-in technique. An advantage, as compared with 
the trapped ion case, is the larger number of photons 
contained by the sidebands for a given n / [23| [30] • 

In summary, we have derived a quantum mechanical 
model of the temperature limit for cooling using radia- 
tion pressure induced dynamical back-action and shown 
that ground state cooling can be achieved as the optical 
cavity linewidth becomes smaller than the mechanical 
frequency, in analogy to atomic sideband cooling. We 
find that the threshold to attain occupancies below unity 
is given by LJ m T > l/y/32. Furthermore, we have shown 
how the spectrum of the optical cavity output could be 
used to measure the final temperatures achieved. Our re- 
sults could apply to other systems exhibiting dynamical 
back-action such as an LC circuit with its capacitance 
modulated by a mechanical oscillator [3]. 
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N+{P) AT_(P )n i + l n f (P) + l C(m + 1)+PA_ 
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